一种从经纬度网格资料获得正方形网格资料的插值方案

本文提出了一种把经纬度网格点上的气象要素值转化到正方形网格点上的多项式插值方案。对插值方案作了稳定性讨论和定量及定性分析。结果表明,这种插值方案方便可行,而且有广泛适用性。     

A study on quality of aquatic environment in Tumen River Area

After the survey of pollution sources, a study on surface water quality assessment and forecast is given by means of grey system method, fuzzy mathematical method and multiple-index method. Based on it, aquatic environment quality features, treatment measures and environmental strategies of the area are proposed. The quality of aquatic environment of 5 rivers in the Tumen River area is studied. The results show that the pollution of surface water is serious; water quality of most rivers is between grade IV and V except the Hunchun River, being higher than grade IV standard; pollution levels of most rivers have been basically controlled except the of Burhatong River, which is deteriorating gradually. Pollutants of the rivers are comparatively regular, mainly are SS, COD, BOD, AR-OH, NH₃-N. The main pollution trades are chemical fibre industry, pulp and paper making industry and mining industry. If the growth rate of gross industrial product is higher than 25 percent under the encouraging-model of regional exploitation, the pollutants’ load will overtake the bearing capacity of aquatic environment. Thus some protection program against pollution must be worked out in order to achieve the harmonious development of economy society and environment. A project from “Studies on Earlier Stage of Regional Development” of the Chinese Academy of Sciences     

Numerical calculations in the orbital determination of an artificial satellite for a long arc

Numerical methods are usually used for the computation of ephemerides with perturbations for the precise orbital determination of an artificial satellite. But their numerical stability will be encountered in a long arc. In this case the use the improved Encke special perturbation methods has been suggested. The results of this paper show that Encke's method does indeed have a certain effectiveness, but cannot yet completely resolve the numerical stability, and the more efficient method is to use the energy integral or its variational relation to control the growth of the along-track error in general numerical calculations so that the aim of stabilization can be achieved.     

雅砻江流域主汛期降水长期预报的一种方法

本文使用天气学方法，从影响雅砻江区域主汛期降水量的长期预报有关因子出发，再运用数理统计方法反复筛选，仅取３个主要因子：西部青藏高原８个站冬春积温，北部冬季极涡面积指数，南部冬季印缅槽平均强度，建立逐步回归预报方程。结果显示：对主汛期６～９月降水量预报拟合准确率为１６／１７即约９４％。回归分析中发现：尽管西太平洋副热带高压对盆地西部降水有很大影响，盛夏且与高原积温有较密切相关，但由放大范围长时段的青藏高原积温，对于雅砻江区域主汛期降水量的影响超过西太平洋副热高压，因而在逐步回归过程中，副高因子仍被剔除．     

武汉区域气象中心数值预报业务并行计算机系统

对武汉区域气象中心并行计算机系统进行了详细地介绍，分析了并行计算机体系结构、网络和存储系统特点；给出了在并行计算机SP上实现数值预报业务并行化的部分结果；对数值预报模式在串、并行编程环境下的结果进行了分析比较。     

A numerical method for differential equations of second order

A new algorithm of order five is presented for the solution of the initial value problem where the system of ordinary differential equations is of second order and does not contain the first derivative.     

Numerical modelling of the impact crater depth-diameter dependence in an acoustically fluidized target

2D numerical modelling of impact cratering has been utilized to quantify an important depth-diameter relationship for different crater morphologies, simple and complex. It is generally accepted that the final crater shape is the result of a gravity-driven collapse of the transient crater, which is formed immediately after the impact. Numerical models allow a quantification of the formation of simple craters, which are bowl-shaped depressions with a lens of rock debris inside, and complex craters, which are characterized by a structural uplift. The computation of the cratering process starts with the first contact of the impactor and the planetary surface and ends with the morphology of the final crater. Using different rheological models for the sub-crater rocks, we quantify the influence on crater mechanics. To explain the formation of complex craters in accordance to the threshold diameter between simple and complex craters, we utilize the Acoustic Fluidization model. We carried out a series of simulations over a broad parameter range with the goal to fit the observed depth/diameter relationships as well as the observed threshold diameters on the Moon, Earth and Venus.     

临猗5.0级地震前后品质因子值变化特征研究

利用西安数字地震遥测台网记录的数字地震资料，采用P波初动半周期残差法求得1998年7月临猗5．0级地震前后不同路径的Q(品质因子)值变化，发现在地震发生前Q值为87～203，震后Q值为67～164，震前震中区附近出现明显的高Q值异常。结果表明，地震前的高Q值异常可以作为地震预测的一种手段。     

Numerical modeling for analyzing thermal surface anomalies induced by underground coal fires

Coal seams burning underneath the surface are recognized all over the world and have drawn increasing public attention in the past years. Frequently, such fires are analyzed by detecting anomalies like increased exhaust gas concentrations and soil temperatures at the surface. A proper analysis presumes the understanding of involved processes, which determine the spatial distribution and dynamic behavior of the anomalies.In this paper, we explain the relevance of mechanical and energy transport processes with respect to the occurrence of temperature anomalies at the surface. Two approaches are presented, aiming to obtain insight into the underground coal fire situation: In-situ temperature mapping and numerical simulation. In 2000 to 2005, annual temperature mapping in the Wuda (Inner Mongolia, PR China) coal fire area showed that most thermal anomalies on the surface are closely related to fractures, where hot exhaust gases from the coal fire are released. Those fractures develop due to rock mechanical failure after volume reduction in the seams. The measured signals at the surface are therefore strongly affected by mechanical processes.More insight into causes and effects of involved energy transport processes is obtained by numerical simulation of the dynamic behavior of coal fires. Simulations show the inter-relation between release and transport of thermal energy in and around underground coal fires. Our simulation results show a time delay between the coal fire propagation and the observed appearance of the surface temperature signal. Additionally, the overall energy flux away from the burning coal seam into the surrounding bedrock is about 30-times higher than the flux through the surface. This is of particular importance for an estimation of the energy released based on surface temperature measurements. Finally, the simulation results also prove that a fire propagation rate estimated from the interpretation of surface anomalies can differ from the actual rate in the seam.     

Three-dimensional discontinuous deformation analysis (3-D DDA) using a new contact resolution algorithm

This paper presents a new contact calculating algorithm for contacts between two polyhedra with planar boundaries in the three-dimensional discontinuous deformation analysis (3-D DDA). In this algorithm, all six type contacts in 3-D (vertex-to-face, vertex-to-edge, vertex-to-vertex, face-to-face, edge-to-edge, and edge-to-face) are simply transformed into the form of point-to-face contacts. The presented algorithm is a simple and efficient method and it can be easily coded into a computer program. In this paper, formulations of normal contact, shear contact and frictional force submatrices based on the new method are derived and the algorithm has been programmed in VC⁺⁺. Examples are provided to demonstrate the new contact rule between two blocks.